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Abstract 

The cosmic string is a useful probe of the early Universe and may give us a clue to physics at 
high energy scales which particle accelerators cannot reach. Although the most promising tool to 
observe it is the cosmic microwave background (CMB), the constraint from gravitational waves is 
becoming so stringent that detecting its signatures in CMB may be impossible. In this paper, we 
construct a scenario that contains cosmic strings observable in the cosmic microwave background 
while evading the constraint imposed by the recent pulsar timing data. We argue that cosmic 
strings with relatively large tension are allowed by diluting loops contributing to the relevant 
frequency range of the gravitational wave background. We also present a particle physics model to 
realize such dilution in the context of chaotic inflation in supergravity, where the phase transition 
occurs during inflation due to the time-dependence of the Hubble induced mass. 
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I. INTRODUCTION 



Ccsndc ..ring^ fl a.e ,ine-,ike topological defects which have been generated in a 

phase transition in the early Universe through the Kibble mechanism [2| . Contrary to 
other topological defects such as monopoles and domain walls, they do not dominate the 
energy density of the Universe, since they evolve according to a scaling law J]- At first, 
they were expected to explain the origin of the primordial density fluctuations that lead 
to structure formation ^|. Although recent observations of the CMB have ruled out the 
possibility that cosmic strings alone are responsible for the primordial density fluctuation, 
they may still contribute to the CMB temperature fluctuation 6|-[9|, which would be a clue 
to exploring physics at high energy scales. 

Thus far, many observations for detecting cosmic strings through the CMB have been 
done by using the Gott-Kaiser-Stebbins effect [10|. The geodesies of photons passing near 
a moving string are perturbed, which affects the temperature fluctuations of the CMB. 
One can thereby constrain the line energy density, or tension, /i. The latest constraint on 
the string tension is Gjj, ^ 10"'' S, Q]*^ with G being the gravitational constant. Future 
observations such as Planck ll[ [l2[ and CMBpol 13 1 are expected to detect the signatures 
of strings on small scales or in larger multipoles of the power spectrum with a tension up to 
G/i ~ 10-8 Is, 9|. 

There is yet another probe to detect or constrain cosmic strings, the stochastic gravita- 
tional wave background (GWB) generated by oscillating string loops 14j-ll9|. When long 
cosmic strings in a Hubble volume intersect with each other, they reconnect and generate 
closed string loops. As a consequence, the distribution of infinite strings obeys the scaling 
rule 



20 



]|-|22|. At the same time, the loops produced in such a way oscillate, emit gravitational 
waves (GWs), and shrink gradually. Since the frequency of GWs is determined by the size 
of cosmic string loops that are generated continuously in the cosmic history, the spectrum of 
GWB, which are the sum of GWs emitted by each loop, are expected to range widely, from 
10^^° Hz to 10^ Hz or even higher frequencies [l^. Although their amplitude is typically 
very small, we may observe their signatures in the GWB through the future gravitational 



*^ See however, for the impossibihty of thermal phase transitions to produce strings relevant to structure 

formation even in the classical big bang cosmology. 
*^ To be precise, the constraint on the string tension depends on the model of string evolution. 
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detectors [19|] such as eLISA|23|], DECIGO [2J], and BBO |25| as well as the ongoing pulsar 



timing arrays (PTAs) 26|, |27|. In PTAs, GWs can be detected through the modulation 



of arrival of pulses from the millisecond pulsars (MSPs) due to the change of the distance 
between the Earth and MSPs caused by the GWB. 

While the PTAs have not detected the signatures of strings thus far, they have obtained 
stringent constraints on the cosmic string tension 26|,|27(|. In particular, the European Pulsar 
Timing Array (EPTA) has recently reported the most stringent constraint, Gfi < 10~^ j^^], 
though it depends on some assumptions regarding string loop formation and GW emission*^ 
Therefore, if we believe this constraint, we cannot hope to detect the signature of cosmic 
strings in the CMB. 

In this paper, we reconsider the formation mechanism of cosmic strings and propose a new 
scenario which does not contradict with the present GW observations while accommodating 
large enough Gfi to be detectable in the CMB. If they are sufficiently dilute, the loop 
formation is reduced and hence the amplitude of GWs is suppressed. This can be realized if 
the phase transition takes place during inflation as proposed in several literatures jsl IsoUsS]. 
While previous proposals assumed nontrivial interactions between the string-forming Higgs 
field and the infiaton field 



3lN33| or gravity 3], [3^, we here point out that the F-term 



infiation in supergravity can naturally realize a phase transition during infiation thanks to 



the Hubble induced mass 
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To be concrete, we construct 
in supergravity 



40|. In 



as studied in Refs. [37|, 
a model that generates cosmic stings during chaotic infiation 
this model, the separation of infinite strings is expanded so large that they cannot make 
loops until long after infiation and hence they enter the scaling regime at a later epoch, for 
example, around the time of matter-radiation equality. Since sources of GWs do not exist 
at early times, we can evade the PTA constraint on the GWB. 

The rest of this paper is organized as follows. In section [Tll we briefiy review the GWs 
emitted by cosmic string loops and see the current limit from observation using MSPs. In 
section IIIIl we propose a new scenario for the generation and evolution of cosmic strings 



Recently, Ref. [28] claims a mild and conservative constraint, G/i < 5 x 10^^, for a certain size of the 
cosmic string loops, which is disfavored by recent numerical simulations [22]. For the size of cosmic string 
loop favored in Ref. [22], Ref. [28] gives a more stringent constraint G/i < 10~^°. Although the size of 
cosmic string loops is still a matter of debate [29] , we here focus on the case where the loop length is large 
as suggested in Ref. [22] and explore a way to evade the PTA constraints. 
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which is consist with current observations and discuss the possibihty of their detection in 
the future. In section llVt we construct a model reahzing such a scenario, and consider how 
to determine the parameters of our model. Section |V] is devoted to discussion. 



II. GRAVITATIONAL WAVES FROM STRING LOOPS AND PTA LIMIT 

Here we review the gravitational wave emission from cosmic string loops. In order to 
estimate the GWB today, first we have to know the size distribution of loops at a given time. 
Without any intersection and reconnection, the energy density of long string network would 
fall in proportion to a~^. Actually, however, these processes frequently occur on cosmic 
time scales to transfer a part of the network's energy to loops which decay by radiating 
gravitational waves if the number density of long strings is large enough. As a result, the 
system relaxes to a scaling solution and the energy density of long strings decreases in 
proportion to d^^ = H^. In terms of the characteristic length of the strings, ^ = 'jdu = 
'-fH^^, where 7 is a numerical coefficient, we can express the energy density of long strings 
as 

p = ^ = j5- (1) 

The coefficient 7 is obtained by numerical simulation. Recent simulation suggests 7 = 0.15 



in the radiation dominant era, and 7 = 0.17 in the matter dominant era (22 1. 

Now let us determine the loop distribution function n{t,i), which represents the number 
density of loops with size £ at time t. At the moment of loop formation [t = tg), its size is 
proportional to the horizon size. We can take the initial length as atg with a a constant. 
The length of a loop becomes smaller with time due to the loss of energy by gravitational 
wave emission. The energy transferred from a loop into gravitational waves per unit time 



is expressed as TGfi^, where F is a constant around 50 2l|]. Thus, the length of a loop at 
time t generated at tg can be written as 

i = atg- TGfiit - tg) . (2) 

In other words, loops with size i at time t were formed at 

_i + TGfd 

^ " « + TG/i ' ^ ^ 
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while those with size i + M sX time t were formed at 

In the flat FRW universe dominated by a perfect fluid with its equation of state given by 

2 

P = wp, the scale factor evolves as a oc t^'^^ , and the Hubble parameter is given by 
2 

H = — — . Taking these facts into account, we obtain the following functional form. 

3{1 + w)t ^ ' ^ 



27(1 + w)^-fHl a{i + TGfit) \ a{t) 



aii + TGfit)^ \a{t) 

The last approximation is valid for loops which were generated in both the radiation- {w = 
1/3) and matter-dominated {w = 0) regimes. 

Now we briefly summarize the formulation to estimate the energy density of GWB. This 
was ori gina lly done by Damour and Vilenkin 15|, and some improvements were proposed 



in Ref. [18|, whose notation we follow. Taking the effect of cusps into account, the Fourier 
mode of amplitude, h, of gravitational radiation emitted by a loop with length i at redshift 
z is given by 

2 

h{f, z, i) = J /r^e[i - euf, z, i)] . (6) 

(1 + Z)3Lpr{z) 

Here ipr{z) represents the comoving distance to the loop defined as 

If' dz' 1 . N 

i^r{z) , (7) 



and 6'ni is defined by = { ^^^2^^ ^ ^ ■ The argument of the Heaviside function O corre- 
sponds to the lowest frequency of loop oscillation. The energy density of GWB is obtained 
by integrating the contribution from various loops, 

^GW(/) = [dz [ dlh\f,Z,i/^^''^^ 



J J ' ^ dzdi 

d'^R(z, tj 

Here — ; — -7^ — is the observable burst rate per length per redshift 

dzdi 

d'R{zJ) _ 1 1 2n{t{z)J)ei{f,zJ) ..... .9^ 



where Hq ^Lpv{z)dz represents the spatial volume corresponding to z z + dz. Its detailed 



derivation is found in Ref. 



18|. 



Next, let us consider the observational constraint on the GWB. Currently the strongest 

limit comes from the Pulsar Timing Array (PTA) experiment, which aims at the detection of 

the GWB using MSPs. (Detailed explanation and eu„ent status ,s found in Ref. fl, TKis 

experiment has sensitivity at / = 10"^ ~ 10^''Hz, which is determined by the duration of 

the observation. At this very low frequency, considering the effect of cusps, one can conclude 

that the GWB spectrum of the energy density, ficw = 7t~~~7^ is proportional to /~^, 

Per a log/ 

coming from loops which are generated during the matter-dominated era. 



According to Ref. [27| the latest limit is 



10-8-4(2a), 



10 



-9.0/ 



(10) 



la). 



(11) 



This has been translated to the upper limit on the line density of cosmic string as [2 

Gfi < 4.0 X 10-^ 

Note that this constraint has been obtained using the method of Damour and Vikenkin [jj 
taking a > TGjj, and incorporating all the relevant modes. *^ On the other hand, analysis 
based on the EPTA data has been done in Ref. 28| incorporating possible variations of 
cosmic string properties. They reported the most conservative constraint on the cosmic 
string tension as 

G/i < 5.3 X 10"^ (12) 
1 value, a ^ 10~^, and assuming the effect 
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18| emphasize that the effect of cusps 



which is obtained in the case of a taking a sma 
of cusps is neglegible. However, recent studies 

which emit higher mode waves, is unignorable. Moreover, recent simulations 
that the typical value of a is much larger than 10~^. If we take the effect of cusps into 



22] indicate 



account and assume a ~ 0.1 as in Ref. 



28 



the constraint of the cosmic string tension turns 
out to be much more severe, < 10"^*^. Here we propose a method by which the EPTA 
constraint can be evaded even in the case where a is not small and when the effect of cusps 
are taken into account. 



The formalism developed in Ref. 18[ would yield a slightly more stringent constraint on the line density. 
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This turns to be possible by considering a different mechanism of string formation. Since 
the dominant contribution to ^Icwif = lyr~^) in the standard scenario comes from the 
loops which were generated during late radiation and early matter-dominated eras, we can 
evade the PTA constraint if the number density of such loops can be reduced. In this case, 
we can explain the PTA constraint without setting Gfi to be small. 



III. A NEW SCENARIO OF STRING FORMATION 

In the standard cosmic string scenario, where strings are formed by spontaneous sym- 



metry breaking after reheating or at the end of hybrid inflation [4l|, |42|] , the string network 
soon relaxes to the scaling solution and it is impossible to reduce the number of small loops 
relevant to the PTA constraint. However, the number of loops can be reduced in a scenario 
where cosmic strings are formed not after inflation but during inflation. This type of scenario 
was also studied by one of us in a different context 

Let us see the detail of this scenario. Since the mean separation of strings becomes larger 
and larger as inflation proceeds, they have to wait for a long time to cross. Strings can 
neither intersect with each other nor form loops until their mean separation falls well below 
the horizon distance. In this scenario, we can remove the sources of the GWB (cosmic string 
loops) at early times, so we can satisfy the observational constraints imposed by GWB with 
a relatively big value of Gfi. 

We calculate ^gw in such a scenario where loops start to form only after the time tsc, 
when strings begin to satisfy the scaling law. Although some loops may have been formed 
before ts^, such loops would be rare and hence we here evaluate the GWB spectrum assuming 
that they are negligible and the GW emission starts at t^c suddenly. 

The calculation can be done as follows. We perform z integral in Eq. (|8]) to z^^, which 
corresponds to the redshift of tgc- The relation between cosmic time t and z is given by 

t = — r — — = —^tiz), (13) 

Ho J. (1 + ^') V^AO + (1 + ^O'^mO + (1 + ^O'^-O ^0 ' ' 

z = ipz{HQt) (the inverse function of (2;)). (14) 

Some modifications are needed in the h integral in Eq. ([8]) . Shortly after tsc, at time ti, the 
minimum loop size is ii = atsc — TGfi{ti — tsc) since very small loops are absent. Besides, 
the B-function means that only the strings with i > ^(^i^^^ = ^2 contribute to the GWB. 



Therefore the h integral runs not from but form max(/ii, /i2), where hi and /12 correspond to 



the amphtude of GWs emitted from the strings with the length ii and £2, respectively. 



?he 
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upper limit of h corresponds to the maximum loop size. Although the authors of Ref. 
insist that one can remove the contribution from rare bursts by introducing the cutoff in 
the h integral, the value of the cutoff is ambiguous. Therefore we do not introduce it and 
include "all bursts". 

Figure [1] shows the density parameter of GWBs from cosmic strings at the frequency 
/ = 1 year~^ as a function of Zsc for various values of a. We can read off the condition 
Zsc — I — 10^ should satisfy for each value of a and Gfi to be consistent with the PTA 
observation. Figure 2 shows the maximum Zsc (or the earliest time for strings to start 
forming loops) in order not to contradict with the PTA data for a = 10~^. We find an 
approximate formula for Zgc, which evades the PTA limit as 

-y \ -0.94 

Therefore, we still have the opportunity to detect features of cosmic string using the future 
CMB observations with small enough ^sc since strings with Gfi ^ 10^^ are allowed if loop 
formation is delayed. 

IV. COSMIC STRINGS AND INFLATION IN SUPERGRAVITY 

We present a new mechanism to realize the above scenario in supergravity. While the 
previous models introduced a nonminimal coupling of the scalar curvature to the string- 
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-|33|, our model does 



forming scalar field js], [s^ or direct coupling with the infiaton 
not require any additional interactions, since it is automatically provided by gravitationally 
suppressed interaction arising from supergravity. 
The F-term potential in supergravity is written as 



K 

V = 



G 



(16) 



where K is the Kahler potential and W is the superpotential. Here we have defined 



K'' = , (17) 

dW 1 dK 

8 





FIG. 1: Energy density of GWB at / = lyear" , I^gW) as a function of Zgc? the redshift 
at the onset of the scahng solution for various values of a. The horizontal line shows the 
upper bound imposed by the PTA observation. 
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The maximum zsc that can evade PTA constraint 




FIG. 2: The maximum that can evade the PTA constraint as a function of G/i for 
a = 10-1. 



and Mg = 1 /VStiG is the reduced Planck mass. During the standard slow-roll inflation, 
the energy density p is dominated by the potential energy, so the Friedmann equation reads 
3H^M^ = p = V. For K < we find 

3M^H^ ^V= {DiW)K'^~\DjWy - -^\W\^ (19) 

in F-term inflation. Since the kinetic terms of the scalar fields are given by —J^*-' (9^ 0^(9^0*, 
the minimal model assumes K = which yields canonical kinetic terms. Here we can 
expand the exponential factor after canonical normalization and diagonalization. 

Therefore, the scalar potential V contains the term 3if^0j0**. This means that the field (pi 
receives an additional Hubble-scale mass, which is called the Hubble-induced mass*^. If the 
Hubble parameter is large enough at the beginning of inflation, these fields fall down to the 
origin quickly and the symmetry is restored during inflation. Moreover, in inflation models in 



*5 



The inflaton should not acquire the Hubble induced mass; otherwise the slow-roll condition cannot be 



satisfied. In order to realize inflation, one must assume shift symmetry 40[ or adopt contrived models 



43| 
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+2 
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Z2 


-1 


-1 












TABLE I: Charge assignments on superfields in the model under the U{l)s symmetry, the R- 
symmetry, and Z2-symmetry. 



which t 
model 



le Hubble parameter changes gradually during inflation, such as the chaotic inflation 



39|, it is possible for the phase transition to take place during inflation. Here we 



adopt the following superpotential. 



W = Wj + Ws, (21) 
Wi = M(f)X, (22) 
Ws = n?s - \s^i), (23) 



and the Kahler potential 



^ = ^(0 + 0*)' + + 1^1' + IV^P + , (24) 

to realize chaotic inflation, where Wj and Ws are superpotentials of the inflation sector 
and the string sector, respectively. Here is the inflaton, X is an additional singlet, -i/^ 
and are the symmetry-breaking Higgs fields that break V{V)s symmetry, s is a singlet 
which destabilizes the Higgs field at the origin, M is the inflaton mass, m is the symmetry 
breaking scale, and A is a coupling constant. Note that we impose i?-symmetry and discrete 
Z2 symmetry in order to suppress all other unwanted interactions such as s0, and so 
on*^. Charge assignments on the fields are shown in Tabled 

This model can realize chaotic inflation naturally thanks to the shift symmetry, — )■ + c 
with c being a real parameter, along the imaginary part of </) which plays the role of the 



*^ Note that the term es^^ in the superpotential cannot be forbidden by these symmetries but it breaks the 
shift symmetry of Since we can expect that the inflaton mass M arises as an order parameter of the 
shift symmetry breaking, the coupling constant e should be suppressed enough, say e ~ A/^ 
Planck units. Thus, the following discussion does not change. 
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inflaton 



40|. Setting = -^(x + ^V^)? the potential for these fields is 



V 



K 



exp 



{DiW)K'^ {DjWy - 



Ml 



\w\ 



Ml 



{x'+\x\'+\s?+m 



Ml 



2 ^ 



(25) 



Note that we do not have to consider all of these terms because the fields Xi 
and ip have very steep potential due to the factor e^'o^ at the field values larger than Mq- 
This means they cannot have field values much bigger than Mq. Since x ^^id s have larger 
masses than the Hubble parameter, they are practically fixed at the origin during infiation. 
Although X has a mass of the order of the infiaton mass, its contribution to the total 
potential energy (or the dynamics of infiation) and the density fiuctuation are negligible due 
to its much smaller field value than the infiaton. Since only the infiaton ip can have a value 
much larger than Mg, this model reduces to a simple chaotic infiation model, V = 

Now we investigate the phase transition during infiation. Diagonizing the mass matrix 
of Ip and Ip, we obtain the following canonically normalized fields: 

\(V^ + V^*), 

(26) 



^2 



V2 



(V' - ^ ) • 



The relevant terms that control the dynamics of the phase transition during infiation are 



V 



2,„2 



Ml V2 



2,„2 



m 



1^1 



ImV + (3if^ - Am^) + 



(27) 



The first term is the potential for the infiaton. The mass term of ipi changes its sign from 
positive to negative during infiation at = Am^/3. As a result, tpi may be destabilized and 
the U{l)s symmetry is spontaneously broken during infiation, while ip2 has always positive 
mass squared. 

Our model contains the parameters A, m, and M. M is related to the primordial power 
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spectrum of curvature fluctuations via 

1 V 1 



247r2 e 



(28) 



k=aH ""^ G 

with A/" ^ 55 being the number of e-folds when observable scales exit the horizon. The 
observed value A^(A;) = 2.4 x 10~^ 4J] determines the inflaton mass to be M = lO^^GeV. 
We assume A is of order of unity and m is much smaller than Mq. 

The line energy density of cosmic strings is related to the vacuum expectation value of 
the string-generating field as *^ 

,2 



/i = 27r(V;/) = 47r— . 



(29) 



Therefore one can determine the value of m? /\ by fixing G/x. 

The crucial issue is whether strings with G/i > 10"*^, whose tension is large enough to 
be detectable by CMB observations |9|, can be produced with appropriate density so that 
they can evade the PTA constraint but are observable using CMB. It is therefore important 
to study when the phase transition takes place to clarify the mean separation and the 
correlation length of the string network. 

In our scenario, the phase transition is triggered by quantum fluctuation during inflation 
rather than thermal fluctuations and we can make use of the results of 




analyze its properties. Calculations based on the stochastic inlfation method [45|, 
that the phase of the string-forming field ipi is fixed when the classical potential force 
surpasses the effect of stochastic quantum fluctuations 35|. This occurs AA/" = a/c/2 e- 



blds after ■j/'i = has become classically unstable with 3if^ = Am^ where c = Am^/M^ 



35 



38 



. The number of e-folds from this epoch to the end of inflation is given by 

Ud = ]:{c-V2c-1), (30) 



which gives the extra dilution which is absent in the case where phase transitions occur after 
or without inflation. The Hubble parameter at this epoch, which we denote with a suffix /, 
is given by 



H 



2c 



1/2 



/ 



M. 



(31) 



*^ To be precise, one should treat the coefficient as a function depending on A, rather than as a constant, 
27r. However within the range of interesting values of A in our model, we can approximate this coefficient 
as 27r[42|. 
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The typical separation between strings at this epoch can be estimated as in [3|, |32|. First, 
choose an arbitrary point Xq with an arbitrary field value of ipilxo). One can always perform 
a gauge transformation to make this real: iPi{xq) = iPir{xq) + iipn^Xo) = ipifj{xQ). Then, 
since a string is a locus of tpu^lx) = tpulx) = 0, the nearby string segment can be found by 
searching for a point (in fact, a line) with ipuilx) = on the surface of ipij{x) = 0. To the 
lowest-order Taylor expansion, the condition for a string to exist within a distance r from 
the point Xq is given by 

r\ei X V^iR(a;o)| > |^iR(a;o)|, (32) 

where ej = ViPii{xq)/\\/iPu{xq)\ is the unit normal vector to the surface ipii = 0. Since 
iPijiIxq) and VtpiR{xo) are statistically independent, almost Gaussian variables, one can 
express the probability for the inequality to be satisfied at a distance r, P(r), in terms of 
the variances (ipfji) = cr^ and {(Vipn)'^) /3 = {{dxipuY) = <Jgi, which will be evaluated later. 
Without loss of generality one can take 2;-axis along ej. Then writing VipiR = {ux,Uy,Uz), 
P{r) is given by 



P{r) = P{^IJiR,Ux,Uy,Uz)d^lJiRd u 



I ^!r Ul + Ul + Ul , 

exp -— — dil}iRd u (33) 



ag>r 



((T2 + a2r2)i/2' 

where P{ipiR, u^, Uy, u^) is the one-point probability distribution function of tpiR and VipiR. 
Thus the typical distance to reach a string, which also gives an estimate of the mean sepa- 
ration of strings, is given by r ^ cr/ag/. 

Next we calculate the variances aX t = tf, which is obtained from the power spectrum of 

^1 



exp 



9 , 2M I (2-V3)2c 



2M2 



H{t 



k) 




(34) 



where 

fi(f^ = ( - - 

2 4 

and tk is the epoch when the comoving mode k exits the horizon. We can evaluate the ratio 



,„.fW-M!_HMV'\ (35, 
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of the variances as 



where 



2 (fk 
(27r)3 



15v^27r(x/ — Xd) exp 



r c 



i7 



1/2 



2c 



2c 



1/2 



^1 



(36) 



(37) 



Now we get the typical separation of strings. In addition to the expression ( l36l) . we also 
numerically evaluate f(c) as a function of c and get f(c) ~ 400 at c = 100. 

After the initial configuration of the string network is determined at during inflation, 
the typical separation simply scales as the scale factor until it falls below the Hubble radius. 
At redshift z, it reads 



d{z) 



re 



Ml 



Or aeq a{z) 



3 X 10^^ 



fic) 



Xf{c) 



exp 



72^-1) 



1 VlO^GeVy 



GeV-\ 



(38) 



where the suffixes end, R, and eq represent the epochs of the end of inflation, reheating and 
radiation-matter equality, respectively. If we assume the network enters the scaling regime 
immediately after the string separation becomes smaller than the characteristic length of a 
scaling network d{z) ^ {(2;), where, 

7 X 7 X 10^1 



-.GeV-\ 



(39) 



iAO -f- fimo(l + Z)^ + ^ro(l + z)- 

then we can get a relation between the model parameres and ^sc by equating (138|) and ( 139|1 
to determine the region of the parameter space where traces of cosmic strings are detectable 
by future CMB observations without conflicting with the PTA constraints. Figure 3 shows 
allowed region of the model parameters where the redshift at the onset of the scaling regime 
has been obtained from an equality 



d{Zsc)/^{Zsc) = 1. 



(40) 



We must admit, however, that this estimate has some uncertainty. Unlike in the standard 
scenario where the string network approaches the scaling solution as the string density de- 
creases, in the present case it is reached as the number of long strings per the Hubble volume 



15 



increases to the scaling value. To our knowledge, such a situation has not been numerically 
simulated, and it is an open issue when the system relaxes to the scaling solution. Fortu- 
nately, however, there exists a finite section of parameter space where all the requirements 
are met, even if we change our criterion of the epoch when the scaling law is achieved. 

Indeed, if we assume the scaling is achieved when an equality d{zsc) / ^ (^sc) = K is satisfied, 
we find that the two vertical lines and the oblique line representing the gravitational-wave 
constraint in Figure 3 shift to the right by 



so that the shape and the size of the allowed region remain intact. Note that the parameter K 
can also accommodate possible variation of the initial string separation, too. Therefore, we 
conclude that there are certainly interestingly large parameter space where signatures of the 
cosmic strings are detectable by the future CMB observation evading the PTA constraints. 

We have shown in Figure 3 the parameter corresponding to Zsc = 1100 and d{z = 0) = 
Hq^ which means there would be only one string inside our present horizon. We find m 
should be about 10^"^GeV, and so use = m/lO^'^GeV. Therefore c becomes lOOAm^^. 
The horizontal axis is related to the time when strings start to scale and produce loops. 
According to the results of the previous section, an upper limit is imposed on Zsc in order 
not to contradict the MSP observation. This upper limit on Zgc yields a lower limit on Xmf^. 
The value of Gfi is proportional to the value of vertical axis. 

V. DISCUSSION 

In this paper, we have studied a new mechanism to evade the PTA constraints on the 
cosmic string tension. We have found that the PTA constraints can be evaded if the loop 
formation starts at a later epoch, say, z ^ 10^, depending on the value of a and the tension 
of the cosmic strings. We have also shown that this scenario can be realized if the phase 
transition and string formation take place during infiation. As a specific example, we have 
shown there is a finite region of parameter space where we can expect to detect the signatures 
of cosmic strings through CMB observations while accommodating the PTA constraints in 
the F-term chaotic infiation model in supergravity j4o| . 




(41) 



2.2 X lO^MogJs: (ate =100) 
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1.00 1.05 1.10 1.15 



FIG. 3: The allowed parameter region of this model is shown. In the blue-painted region, the 
network of strings starts to produce loops so early that the gravitational waves emitted by them 
contradict observations. Though determining this region requires more statistical analysis of 
bursts and a more realistic treatment of a, we emphasize that some region remains allowed. 
The most interesting region, which represents the case strings in which can not only evade 
the PTA constraint but also contribute to the anisotropy of the CMB as much as ordinary 
strings do. The white region corresponds to the case where strings are observable through the 
CMB but at densities smaller than in the conventional scenario. As a recent constraint on 



Ql- According to ^ 



G/i, we used the value obtained in [7|]. According to 
observable. 

As seen in Figure 1, t 



strings with G/i < 7 x 10 ^ are not 



;he value of a is very important to determine the constraint on 
Zsc- Recent simulation (22| indicates that loops are formed with various sizes ranging from 
a ~ 10~^ to a < 10~^ in the standard scenario. In our scenario, however, the mean 
separation of cosmic strings is much larger than the horizon initially, so it is expected 
that larger loops are produced more frequently than in the ordinary scenario. In addition, 
inflation also erases small-scale fluctuations on long strings. Hence strings do not develop 
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small-scale structure and the correlation length remains of the order of the Hubble radius for 
some time, even after Zgc- As a result, we expect a to take a relatively large value, a ~ 10~^, 
at the onset of the scaling regime. 

Though we assume that the cosmic string network enters the scaling regime when the 
mean separation of cosmic strings becomes comparable to the characteristic length of the 
scaling network, the time of the onset of the scaling regime is still ambiguous. However, 
we can expect that the earliest candidate for the onset of the scaling regime is when the 
smallest scale of the cosmic strings, (a{z) /af)Hj^ , becomes comparable to the characteristc 
length of the scaling network ^{z). We confirmed that, even in this case, there is a finite 
but interesting parameter space. Therefore, we conclude that there is certainly a region of 
parameter space where we can expect to detect the signatures of cosmic strings by future 
CMB observations. 

In our scenario, the effect of the CMB may need to be reconsidered. If inflation did not 
last so long after strings were formed, stings start to satisfy the scaling rule at an earlier 
epoch. As a consequence, their separation would already become smaller than the horizon 
at the recombination. Then there would be as many strings as in the standard scenario 
near the last scattering surface of the CMB photons. The effect on the CMB would be the 
same as that has already been analyzed in the literature 47|]. On the other hand, if inflation 
lasted very long after strings had been formed, they could not start to scale until very late 
time. Therefore, their separation would be larger than the horizon at the recombination. In 
this case, their effect on the CMB would become smaller. For a quantitative understanding, 
numerical simulations of the string network with appropriate initial conditions are desired. 
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